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ABSTRACT 

An exploration of the relationship between bolometric luminosity and outflow velocity for 
two classes of X-ray outflows in a large sample of active galactic nuclei has been performed. 
We find that line radiation pressure could be one physical mechanism that might accelerate 
the gas we observe in warm absorber, v ~ 100-1000 km s _1 , and on comparable but less 
stringent grounds the ultrafast outflows, v ~ 0.03-0. 3c. If comparable with the escape velocity 
of the system, the first is naturally located at distances of the dusty torus, ^lpc, and the 
second at subparsec scales, MX01 pc, in accordance with large set of observational evidence 
existing in the literature. The presentation of this relationship might give us key clues for 
our understanding of the different physical mechanisms acting in the centre of galaxies, the 
feedback process and its impact on the evolution of the host galaxy. 
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1 INTRODUCTION 

Mildly relativistic and non-relativistic absorption troughs are ob- 
served in the X-ray spectra of active galactic nuclei (AGNs). A 
qualitative separation is usually done between classical v — 100- 
1000km s -1 warm absorbers, which we refer to (throughout this 
Letter) as low-velocity outflows (Blustin et al. 2005; McKeman, 
Yaqoob & Reynolds 2007), and ultrafast outflows (UFOs) v > 
10 000 km s _1 , which we refer to as high-velocity outflows (Tombesi 
et al. 20 10a, b). 

On one hand, the extracted velocity from the classical v ~ 100- 
1000 km s” 1 warm absorber is mainly obtained using the Fe M- 
shell 2p-3d unresolved transition array (Behar, Sako & Kahn 2001; 
Ramirez et al. 2008) and O vn or O vui resonance lines (e.g. George 
et al. 1998; the hallmark of the classical warm absorber). The spatial 
location of this absorbing material is uncertain. Models of X-ray 
absorbers in AGN place them at a wide range of distances from 
the central source. Specifically, they are suggested to be winds 
originating from the accretion disc (Murray et al. 1995; Elvis 2000), 
located at the dusty (M pc) torus scales (Krolik & Kriss 2001) or 
even beyond the narrow-line region (e.g. Ogle et al. 2000). 

On the other hand, blueshifted Fe K absorption lines have been 
detected in recent years at E > 7 keV in the X-ray spectra of sev- 
eral radio-quiet AGNs (e.g. Chartas et al. 2002; Chartas, Brandt & 
Gallagher 2003; Pounds et al. 2003; Dadina et al. 2005; Markowitz, 
Reeves & Braito 2006; Braito et al. 2007; Ramirez 2008; Cappi et al. 
2009; Reeves et al. 2009; Tombesi et al. 20 10a, b). They are usually 
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interpreted as due to resonant absorption from Fe xxv-xxvi associ- 
ated with a zone of circumnuclear gas photoionized by the central 
X-ray source, with ionization parameter log £ ~ 3-6 erg s _1 cm and 
column density Nu ~ 10 22 ~10 24 cm“ 2 . The energies of these ab- 
sorption lines are systematically blueshifted and the corresponding 
velocities can reach up to mildly relativistic values of ~0.03c-0.3c. 
These findings are important because they suggest the presence of 
previously unknown massive and highly ionized absorbing material 
outflowing from their nuclei, possibly connected with accretion disc 
winds/ejecta (e.g. King & Pounds 2003; Proga & Kallman 2004; 
Sim et al. 2008, 2010a; Ohsuga et al. 2009; King 2010b). 

Several acceleration mechanisms have been proposed to explain 
these outflows: (1) thermally accelerated winds (Krolik & Kriss 
2001); (2) radiation pressure through Thomson scattering and mag- 
netic forces (MHD; Ohsuga et al. 2009) and (3) radiation pressure 
due to the absorption of spectral lines (e.g. Proga & Kallman 2004; 
Ramirez 2008, 2011; Schurch, Done & Proga 2009; Sim et al. 
2010b). Although the first one can explain the velocities we ob- 
serve in the low-velocity outflows, it can be excluded because it 
cannot explain the ~0.1c-0.2c we observe in UFOs (Tombesi et al. 
2010a,b). Ohsuga et al. (2009) seem to reproduce the velocities ob- 
served in low - and high-velocity outflows. On the other hand, Arav, 
Li & Begelman (1994), Ramirez (2008), Saez, Chartas 8c Brandt 
(2009) and Chartas et al. (2009) invoke radiation pressure due to 
lines to explain the ~0.2c outflow detected in a good signal-to-noise 
ratio X-ray spectrum of a high-z quasar (the broad absorption line, 
APM 08279+5255), and they reproduce, as part of the procedure, 
the Fe xxv-xxvi lines detected at E > 7 keV. Here we focus on this 
kind of approach, since it allows us to explore, very efficiently, a 
wide range of physical parameters of the system. 
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The goal of this Letter is to place a proof of idea for a sys- 
tematic study about the operating acceleration mechanisms in both 
X-ray low - and high-velocity outflows, using an anisotropic radia- 
tive pressure framework (e.g. Proga, Stone & Kallman 2000; Proga 
& Kallman 2004; Liu & Zhang 201 1), beginning with line radiation 
pressure. 

We present the observables from which we build the model 
in Section 2. The details of the proposed model are presented in 
Section 3. The results and the discussion are in Section 4. We sum- 
marize in Section 5. 


2 OBSERVABLES 

In this section we present the observables of the two types of out- 
flows, since they are the initial motivation for the proposed model. 


2.1 The low-velocity outflows 

When describing the physical conditions of warm absorbers, it is 
common to use the definition of ionization parameter f 
(Tarter, I\icker & Salpeter 1969), where F ion is the total ionizing 
flux (F ion = L ion /47tr 1 2 ) and ah is the gas density. The source 
spectrum is described by the spectral (specific, energy-dependent 
e) luminosity L € = L i0l / o where Li 0n is the integrated luminosity 
from 1 to 1000 Ryd, and J* mRyi f € de = 1. 

So we describe the ^1000 km s _1 warm absorber outflows as ab- 
sorbing material around a supermassive black hole (SMBH) with 
mass M B h ~ 3 x 10 7 (Peterson et al. 2004; Blustin et al. 2005), col- 
umn density of the absorbing material iV H ~ 10 20 — 10 22 cm -2 , flow- 
ing outwards at velocities v ~ 100-2000 km s _l (e.g. Kaspi et al. 
2002; Krongold et al. 2003, 2005; Ramirez, Bautista & Kallman 
2005), at medium ionization states log £ ~ 0-3 erg s” 1 cm (Blustin 
et al. 2005). 

When computing the energetics, they find mass-loss rates of 1 
M out ^ 0.6 Mq yr 1 , ratios of Af out to accretion rates 2 M out /Macc w 
5 and kinetic luminosity 3 L EK = of the orders of 10 38 - 

10 41 ergs^ 1 , representing less than 1 per cent of the bolometric 
luminosity (Blustin et al. 2005). The main conclusion from these 
estimations is that these outflows contributes little to the energy 
injected in to the host galaxy. But the amount of matter processed 
over the AGN lifetime can be significant (also in accordance with 
Krongold et al. 2007, for instance). 


2.2 The high-velocity outflows 

The characteristics of the UFOs with v > 10 000 km s -1 (>0,03c) 
can be derived from the blueshifted Fe xxv-xxvi absorption lines 
detected by Tombesi et al. (2010b) in a complete sample of local 
Seyfert galaxies. Such features are detected in ~40-60 per cent 
of the sources, which suggests a covering fraction of C ~ 0.5. 
Tombesi, Cappi & Reeves (2011) also performed a photoionization 
modelling of these lines. They derived the distribution of the outflow 
velocities, which ranges from ~0.03c up to mildly relativistic values 
of ~0.3c, with a peak and mean value at ~0. 14c. As expected, these 


absorbers are highly ionized, with log £ ~ 3-6 erg s -1 cm, and have 
large column densities, in the range iV H ^ 10 22 — 10 24 cm -2 . 

The SMBH masses of the Seyferts in the Tombesi et al. (2010b) 
sample have a mean value of M B h = 5.3 x 10 7 Mq (Marchesini, 
Celotti & Ferrarese 2004; Peterson et al. 2004). 

When computing the energetics of these outflows, we see 
that these are more massive than the low-velocity ones, M out ~ 
0.1-1 Mqjt - 1 ~ Mace, and also much more powerful, with a 
mechanical power of ~10 43 -10 45 erg s -1 (e ; g. Pounds et al. 2003; 
Markowitz et al. 2006; Braito et al. 2007 ; Cappi et al. 2009; Tombesi 
et al. 2010b). The latter value is ~5-10 per cent of the bolometric 
luminosity. Therefore, the high-velocity outflows may potentially 
play an important role on the expected cosmological feedback from 
AGNs (e.g. King 2010a, b). 


3 THE MODEL 


We build our model based on observational evidence, basically 
from two sources: Blustin et al. (2005) for the classical v ~ 100- 
1000 km s -1 warm absorbers, referred here as low-velocity outflows, 
and Tombesi et al. (2010b) for the UFOs v > 10 000 km s -1 , referred 
here as high-velocity outflows. 

As we see, together the two classes of outflows cover a wide 
range in velocity and black hole masses. We seek a physical model 
which provides the context to explain (to first approximation) part 
of the observational evidence we have up to now. 

In order to do that, and to gain some insight into the relation- 
ship between outflow velocity (u out ) and luminosity, we invoke the 
velocity profile as a function of radius, given by hydrodynamical 
calculations (Proga et al. 2000). The mathematical shape of the 
relation is given by assuming an outflow accelerated by radiation 
pressure from a central source with a bolometric luminosity 
and a mass of M B h- 


iw[*'] - 


2GM B h[«] (Vf—lfl 
V ^Edd [l] 




( 1 ) 


where Lem is the Eddington luminosity, T f is the force multiplier, 
where the acceleration due to the absorption of discrete lines is 
encapsulated (Laor & Brandt 2002), R ^ is the radius at which the 
wind is launched from the disc, and R is. the distance of the accel- 
erated portion of the outflow from the central source. The index i 
runs simultaneously over the two distributions: velocity/BH mass. 
Assuming that we observe the gas when it has reached the terminal 
velocity of the wind, i.e. u out at R = oo, we can write 
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Now we have the basic ingredient of the model, and we are ready 
to build our two classes of outflows. 

Set 1: high velocity. This synthetic sample is composed from 1000 
black holes with masses normally distributed with mean /x m = 2.6 x 
10 7 Mq and standard deviation o m = 1.3 x 10 7 Mq . 4 Also we use 
a normal distribution for the outflow velocity of the absorbing gas 
around the BH with = 57 000 km s^ 1 (this is the best-fitting value 

for set 1) and cr v = 15 000km s -1 . 


1 Mean of the Mmi reported by Blustin etal. (2005) in their table 4, excluding 
Ark 564 (outlier M 0 ut = 23 Mq yr~ l ). 

2 Mean of the M 01 u/Af acc reported by Blustin et al. (2005) in their table 4, 
excluding Ark 564 (outlier Afout/Micc = 550). 

3 Excluding PG 0844+349 and PG 1211+143. 


4 We take the middle point (jt m = 2.6 x 10 7 Mq) between the mean ex- 
tracted from Tombesi et al. (2011) (/u. m — 5.3 x 10 7 Mq), excluding the 
mass of Mrk 205 (outlier M bh = 44 x 10 7 Mq), and the mean extracted 
from table 5 of Blustin et al. (2005) (^ m = 2.7 x 10 7 Mq), excluding the 
mass of IRAS 13349+2438 (outlier Mbh = 80 x 10 7 Mq). Also we take 
the larger value of the observed dispersion cr ra (obs) = 1.3 x 10 7 Mq. 


© 2011 The Authors, MNRAS 419, L64-L68 

Monthly Notices of the Royal Astronomical Society © 2011 RAS 



L66 J. M. Ramirez and F. Tombesi 



Figure 1. Velocity of the outflow versus the bolometric luminosity needed 
to accelerate the wind. Filled circles are the velocities observed if we see 
the wind along the flow. Open squares are including the effects of random 
lines of sight (see text for discussion). 


Set 2: lov velocity . In this case we use 1000 black holes with masses 
normally distributed with the same mean, /x m = 2.6 x 10 7 Mq, and 
standard deviation, cr m = 1.3 x 10 7 Mq, as before but we use a 
normal distribution for the outflow velocity of the absorbing gas 
around the BH with jx v = 1800 km s _l (this is the best-fitting value 
for set 2) and c v = 600 km s _l . 

The best-fitting values of /jl v ( 1) and fi v (2) [as well as Ff(l) = 250 
and r f (2) = 40; see next section] are estimated by comparing a grid 
of models computed using 10000 < /x*(l) < 70000 (kms -1 ) and 
100 < (2) < 2000 (km s _l ), with the corresponding set of obser- 

vation: set 3 versus Tombesi et al. (2011) and set 4 versus Blustin 
et al. (2005). We performed the comparison using a Kolmogorov- 
Smimov test. We take the /x„(l) and jx„(2) [as well as Tf(l) = 250 
and Tf (2) = 40; see next section] which give maximum p- values of 
the test. The final p-value for the comparison set 3 versus Tombesi 
et al. (20 1 1) gives ^0.8 (D = 0. 16) and for the set 4 versus Blustin 
et al. (2005) M).07 (D = 0.38). 

The launching radius, R m > for each set of simulations is based 
on the results of Blustin et al. (2005) for set 2, i.e. /? in (set2) = 1 pc 
(orders of magnitude value). For set 1, we set /^(setl) = 0.01 pc, 
based on Tombesi et al. (20 10a, b). 


4 RESULTS AND DISCUSSION 

In this section we discuss the results of our simulated samples and 
summarize their physical context and implications. 

Fig. 1 shows the theoretical relationship between outflow velocity 
and bolometric luminosity. 

The curvature we see in v out of set 2 (and also set 1 , filled circles, 
our line of sight is along the flow; see below) is due to the quadratic 
dependence of the luminosity with v out given by equation (2). It is 
worth noting that if we assume a force multiplier of T f = 40 (this is 
the best-fitting value for set 2 and set 4; see Section 3), we are able 
to reproduce the range of luminosities seen in fig. 4 of Blustin et al. 
(2005). 

On the other hand, this model is not reproducing well the dis- 
persion we observe in fig. 4 of Blustin et al. (2005), which could 
be explained by the fact that equation (2) will give radially acceler- 
ated flows, or in other words, that we are observing all the objects 
radially along the flow. To include the effect of observing random 
lines of sight transversely through different sections of the flow, we 
connect [he observed velocity v 0 b S with the intrinsic radial veloc- 
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Figure 2. Theoretical versus observational outflow velocity against lumi- 
nosity. Ellipses (la, 2a and 3a contours of the model) are theoretical 
calculations coining from two samples of outflows - left-hand panel: high 
velocity (set 3, i.e. including the effects of random lines of sight) and right- 
hand panel: low velocity (set 4, i.e. including the effects of random lines of 
sight). Open squares are observational points from the XMM-Newton radio- 
quite sample of Tombesi et al. (2011), for the UFOs. The open diamonds 
are points compiled by Blustin et al. (2005). 


ity, i.e. t> 0 bs = v O utCOS0, where 0 is the angle between the outflow 
direction and the line of sight. Then using the same luminosities 
produced by set 2, but plotting (open squares) against 1000 u 0 bs(s) 
randomly computed using a random generator of numbers (with 
ti/12 < 0 < tt/ 2), 5 and equation (2) with u 0 bs = v ovt cos0, we 
produce the subsample set 4. It is easy to see that set 4 resembles 
better the plot shown in fig. 4 of Blustin et al. (2005), taking into 
account the possible incompleteness of the sample. 

Doing the same for the UFOs, we use the luminosities produced 
by set 1 (Tf = 250, this is the best-fitting value for set 1 and set 
3; see section 3) but plotting (open squares) against 1000 t> 0 bs(s) 
randomly computed using a random generator of numbers (with 
7T/12 < 0 < tt/2), and equation (2) with u 0 bs = iwcosfl, to 
produce the subsample set 3. In this case, set 3 better resembles the 
relationship between v^t and luminosity (see Fig. 2 ). 

In Fig. 2, we place both classes of outflows together along with 
observational points taken from two samples of objects: 14 points 
out of the 23 objects reported by Blustin et al. (2005) (the others 
either did not report outflow velocity or were unknown), in the right- 
hand panel, and 19 points from the sample of 19 objects where UFOs 
have been detected by Tombesi et al. (201 1), in the left-hand panel. 

There are several interesting facts about the plot: (1) the obser- 
vational points cover «3.8 orders of magnitude in velocity; (2) 
they cover ~3 orders of magnitude in bolometric luminosity (from 
-lO 43 to 10 46 ergs -1 ) and (3) that our proposed model is able to 
reproduce most of the low- velocity points (1 1/14) (two of the points 
might fall in the high-velocity set instead) and, less stringently, half 
of the points for the high-velocity set (11/19), using one physi- 
cal acceleration mechanism and three free parameters: mass of the 
BH (A/bh)» outflow velocity (iw) and force multiplier (T f ). On 
the other hand, there are two sectors where the deviations between 
model and theory are large: (1) the high velocity/low luminosity 
and (2) the low velocity/high luminosity, both requiring a closer 
inspection to the completeness of the samples and/or the addition of 
other acceleration mechanisms, like magnetic thrust, for instance. 
But, this may be the topic of future work. 


5 The lower limit assumes that the torus covers ~30°, so we are able to 
observe the outflow only from 6 > tc/12. 
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4.1 Anisotropic radiation pressure picture 

We propose here that the low- and high-velocity outflows can be 
accelerated by the same physical mechanism (as is suggested by 
Fig. 2), i.e. radiation pressure, and that the differences depend 
on the values of the three fundamental parameters: mass of the 
BH (MbhX the object luminosity (L b oi) and force multiplier (r f ). 
The values ofTf we have used are of the orders of those found in 
observational (e.g. Laor & Brandt 2002) and theoretical (e.g. Saez 
et al. 201 1) works. However, detailed photoionization computations 
are required to verify if the opacity of the gas under the physical 
conditions presented here can overcome the overionization problem 
(Proga et al. 2000), and are left for future work. 

The anisotropic property of the radiation is basically demanded 
by the existence of obscured and unobscured (Type 1 and Type 2) 
AGNs (Antonucci 1993), and it explains the decreasing of Type 2 
AGNs as a function of the X-ray luminosity (Hasinger 2008). It 
is also intrinsically linked to the existence of the dusty torus 6 and 
it gives the natural frame to locate our two classes of outflows. If 
comparable with the escape velocity of the system (a black hole of 
Mbh = 5.3 x 10 7 Mq), fast outflow, v oux ~ 0.10c, locates escape 
radius at r ~ (20-100)r s (Schwarzschild radius, r s = 2G ff BH )• This 
is very close to the SMBH, and the origin is likely the accretion 
disc. 

Again, using 500-1500 km s -1 as escape velocity from a 2.6 x 
10 7 Mq BH locates the gas at ~0.1— 1 pc, well beyond the event 
horizon of the BH, and the broad-line region (Laor & Draine 1993) 
as well, where the dusty torus is thought to be located. 

5 SUMMARY 

For the first time, an exploration on the relationship between bolo- 
metric luminosity and outflow velocity for two classes of X-ray 
outflow; in a large sample of AGNs has been performed. We find 
that (1) line radiation pressure is an efficient 7 mechanism to ac- 
celerate the low-velocity (500-2000 km s _1 ) gas we observe in the 
classical MOOOkms -1 warm absorber, (2) it might also become 
efficient to accelerate the high-velocity (0.05-0,3c) gas we observe 
in the UFOs and (3) they both might be placed in the same context 
of anisotropic radiation pressure (e.g. Proga et al. 2000; Proga & 
Kallman 2004; Liu & Zhang 2011). 

However, there are many open questions, which will require close 
investigation and detailed modelling. In the first place, the fact 
that we are assigning one type of outflow to different portions 
of the parameter space (M BH , Lboi and T f ) does not preclude the 
existence of both in the same object. Also, careful studies of the 
connection between these outflows observed in X-rays and outflows 
seen in other bands of the electromagnetic spectrum, ultraviolet, 
infrared or optical (multiwavelength studies) are important points 
to be addressed. Finally, the close inspection of the main parameters 
of the system (A/bh, Worn and T f ) with cosmological parameters (like 
redshift) is of high relevance, and might be the topic of future works. 
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